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1.  INTRODUCTION 


The  adiabatic  auroral  arc  model  [Chiu  and  Schulz,  1978;  Chiu  and  Cornwall, 
1980;  Chiu  et  al.,  1981;  Lyons,  1980,  1981]  is  a  reasonable  if  incomplete 
description  of  a  time-independent  auroral  arc  driven  by  a  large-scale  parallel 
potential  drop  (as  seen  experimentally;  see,  e.g.,  Reiff  et  al.  [1988]).  Even 
the  static  properties  of  the  model  are  difficult  to  investigate  because  the 
defining  equations  are  nonlinear,  but  recently  Cornwall  [1988]  has  found  some 
exact  two-dimensional  (latitude  and  longitude;  height  dependence  suppressed  by 
height  integration)  stationary  solutions. 

The  next  step  in  investigating  the  model  is  to  look  at  the  low-frequency 
dynamics,  which  is  expected  to  show  unstable  Kelvin-Helmholtz  [e.g.,  Hallinan 
and  Davis,  1970]  modes  and  possibly  other  interchange  modes  as  well.  Several 
physically  different  time  scales  are  contained  in  the  adiabatic  model,  and 
they  all  turn  out  to  be  in  the  range  of  tens  of  seconds: 


Ionization  growth  time  T^ 
E  x  B  drift  transport  time  Tp 
Dissociative  recombination  time  TR 


eN  _ 

roa*  “  30  3ec 

(1.1) 

BA2 

— —  s  50  sec 
cA$ 

(1.2) 

(2aN) -1  *  50  sec 

(1.3) 

In  these  equations,  c,  e,  and  B  are  the  speed  of  light,  the  electronic  charge, 
and  the  Earth's  magnetic  field,  respectively,  and  the  other  quantities  and 
their  nominal  values  are: 

12  -2 

Height -integrated  ionospheric  plasma  density  N  *  10  cm 
Number  of  ion  pairs  per  incident  auroral  primary  r  *  30 

Current -potential  parameter  Q  *  0.1 -0.2  cm  1  sec  *  (see  (1.4)  below) 

Nominal  parallel  potential  drop  A0  a  1 . 5  kv 

14  -2  -1 

Height -integrated  dissociative  recombination  rate  a  2  10  cm  sec 
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Invert ed-V  scale  length  A 


*  40  km 


(see  (2.7)  below) 


Many  of  these  parameters  can  deviate  substantially  from  their  nominal  values, 
especially  as  the  parallel  potential  drop  changes,  and  so  one  or  another  of 
the  three  time  scales  in  (1.1,2, 3)  may  dominate  (i.e.,  be  smallest).  For 
example,  a  low-energy  aurora  may  not  penetrate  to  the  heights  where 
dissociative  recombination  is  important,  so  a  will  be  much  smaller  then 
nominal.  Or  just  after  an  aurora  turns  on  in  the  midnight  ionosphere,  N  may 
be  an  order  of  magnitude  smaller  than  nominal,  and  again  recombination  is 
unimportant.  In  writing  the  nominal  parameter  values  we  have  not  attempted  to 
construct  an  internally  consistent  set;  later  in  applications  we  will  be  more 
precise. 

There  are,  of  course,  a  host  of  other  time  and  space  scales  that  one  might 
invoke  in  auroral  dynamics  associated  with  local  wave  modes  :  (drift,  ion 
cyclotron,  lower  hybrid,  .  .  . )  or  with  magnetospheric  phenomena  (wave,  particle 
quarter-bounce  times,  polarization  currents.  Pc  1-5  waves,  reconnection  time 
scales) .  Also  there  are  correspondingly  marty  studies  of  phenomena  associated 
with  these  scales,  as  well  as  the  ones  that  concern  us;  for  example,  Basu  et 
al.  [1988],  Gorney  [1989],  Weimer  et  al.  [1985,  1987]  study  the  spatial  scales 
of  auroras  from  0.1  to  100  km  in  the  ionosphere  or  magnetosphere.  Theories  to 
explain  the  resulting  power  spectra  are  reviewed  in  Kintner  and  Seyler  [1985]; 
Keskinen  and  Ossakow  [1983];  Fejer  [1989];  and  Huba  [1989].  In  most  cases 
theoretical  analysis  does  not  go  beyond  linear  instabilities,  and  is  then 
supplemented  with  large  computer  simulations,  such  as  that  of  Mitchell  et  al. 
[1985]  for  high-latitude  F-layer  instabilities  and  of  Keskinen  et  al.  [1988] 
for  Kelvin-Helmholtz  (KH)  instability  in  the  high-latitude  ionosphere.  Linear 
analysis  is  complicated  by  the  need  to  account  for  the  destabilizing  effects 
of  auroral  gradients,  whether  on  large  scales  for  such  modes  as  Kelvin- 
Helmholtz,  or  on  small  scales  (e.g.,  Ganguli  et  al.  [1985];  Basu  and  Coppi 
[1989]).  The  particular  mechanisms  we  emphasize,  namely  a  parallel  current 
driven  by  a  parallel  potential  drop  and  the  consequent  ionization  of  the 
auroral  ionosphere,  also  have  consequences  on  spatial  scales  larger  than  the 
typical  transverse  width  of  an  inverted-V,  such  as  current  oscillations  on 
scales  of  hundreds  of  km  [Weimer  et  al.,  1988]  and  modification  of  large- 
scale  convection  fields  (e.g.,  Blomberg  and  Marklund  [1988]).  Finally,  these 
mechanisms  can  react  on  the  driving  magnetospheric  source,  causing  feedback 
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effects  which  we  do  not  attempt  to  analyze  here  [Lotko  et  al.,  1987], 

On  the  one  hand  we  would  like  to  incorporate  all  these  effects  into  a 
comprehensive  model  of  time-dependent  auroral  dynamics,  but  so  far  that  is  out 
of  the  question.  On  the  other  hand,  there  is  as  yet  little  theoretical  work 
and  no  simulations  which  include  the  specific  effects  mentioned  above  which 
are  the  ».  ilimarks  of  the  adiabatic  auroral  model.  The  first  is  an  auroral 
parallel  current  J  ^  (positive  upward)  whose  strength  is  linear  in  the 
difference  in  potential  between  the  ionosphere  (0)  and  the  equator  (0  ) : 

J„  =  Q ( 0  0 E )  •  (1-4) 

The  second  is  an  auroral  ionization  source,  of  strength  Tj||/e,  where  T  is  the 
number  of  electron-ion  pairs  produced  by  a  single  auroral  primary.  In 
principle,  T  itself  depends  linearly  on  0-0  (since  T  is  the  ratio  of  the 

Cj 

auroral  primary  energy  to  the  energy  needed  to  form  a  pair) ,  but  for  the  most 
part  in  this  report  we  will  take  T  to  be  a  given  constant;  this  does  not  affect 
any  of  our  qualitative  conclusions  (cf .  Cornwall  [1988]  who  studied  both  cases 
for  T  for  the  static  auroral  arc  model)  .  The  ionization  sources  we  use  in 
this  report,  and  their  significant  effects  on  auroral-ionospheric  densities 
that  we  find,  are  consistent  with  experimental  studies  of  spatially  structural 
auroral  ionization  produced  by  precipitating  electrons  [Labelle  et  al.,  1989] 

The  purpose  of  the  present  report  is  to  study,  in  an  analytical  way,  the 

influence  of  the  ionization  source  represented  by  when  recombination  is 

insufficient  to  bring  about  equilibrium;  that  is,  when  T  of  (1.3)  is  long 

R 

compared  to  Tr  of  (1.1).  This  can  happen,  in  particular,  for  low-energy 
auroras  which  deposit  their  energy  at  altitudes  where  dissociative 
recombination  is  small.  One  might  then  expect  that  transport  processes  would 
take  over,  yielding  an  equilibrium  where  ionization  is  transported  to  the 
return-current  region  a3  fast  as  it  is  produced  in  the  aurora.  But  we  show 
that  this  is  generally  impossible;  there  are  always  regions  where  horizontal 
transport  vanishes  or  is  very  small.  There  are  then  at  least  three 
possibilities  to  control  the  growth  of  ionization  in  such  a  region:  (1)  the 
region  of  small  horizontal  transport  moves  elsewhere;  (2)  vertical  downward 
transport  to  altitudes  where  recombination  is  effective  takes  place;  (3)  the 
ionospheric  potential  0  decreases,  thus  lowering  the  ionization  source 
strength  (see  (1.4)). 
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In  this  report,  we  concentrate  on  possibilities  (1)  and  (3)  above,  and  we 
will  find  that  they  indeed  are  important.  As  for  possibility  (2),  vertical 
transport,  the  projection  of  neutral  winds  along  B  in  the  auroral  zone  is  too 
small  to  be  important,  and  the  vertical  component  of  such  winds  in  the  aurora 
is  upward  (because  of  auroral  heating) .  Upward  winds  act  in  the  wrong 
direction  to  bring  excess  ionization  downward,  but  they  can  raise  molecular 
ions  to  greater  altitudes,  chus  promoting  dissociative  recombination  there. 
But  the  upward  velocity  is  only  a  few  m/sec,  too  small  to  matter.  It  is  true 
that  the  upward  vertical  winds  are  deflected  into  horizontal  flow  which  would 
carry  ionization  out  of  the  aurora,  in  a  horizontal  flow  pattern  which  appears 
to  have  a  strong  divergence,  but  the  wind  speeds  in  these  cells  is  only  tens  of 
m/sec,  too  small  to  be  important. 

We  will  analyze  the  problem  of  ionization-driven  auroral  dynamics  in  two 
ways.  The  first  is  traditional  local  perturbation  theory,  which  yields  a  long- 
wavelength  growth  rate  with  a  term  proportional  to  the  ionization  rate  T. 
Other  terms  enter  at  shorter  wavelengths  which  we  identify  with  corresponding 
terms  in  the  usual  E  x  B  -gradient -drift  instability;  one  expects  that  in 
general  these  transport  terms  tend  to  damp  the  instability,  but  this  need  not 
always  be  the  case.  Finally,  there  are  terms  which  represent  the  direct 
influence  of  the  parallel  current  of  (1.4).  Their  effect  is  not  comparable 
to  parallel  currents  in  the  current-driven  E  x  B  -gradient-drift  instability, 
where  a  parallel  current  perturbation  involves  a  density  perturbation. 
Evidently  from  (1.4)  a  perturbation  in  Jjj  is  a  perturbation  in  <J>,  which  has 
quite  different  consequences.  It  is  more  appropriate  to  describe  our 
instability  as  an  ionization-driven  E  x  B  -gradient -drift  instability. 

Second,  we  are  able  to  find  3ome  locally  exact  solutions  to  the  full 
nonlinear  equations  of  the  model.  We  call  them  "locally  exact"  because  while 
they  are  indeed  exact  solutions,  the  assumed  spatial  dependence  (namely, 
quadratic)  of  the  height-integrated  density  N  and  the  potential  0  is  not 
sufficiently  general  to  admit  extension  to  the  entire  auroral  region  out  to 
the  return-current  boundary.  (There  are  one  or  two  relatively  uninteresting 
globally  exact  time -dependent  solutions,  which  we  mention  in  Section  5.) 
These  locally  exact  solutions  show  the  phenomena  mentioned  above,  in  which  the 
regions  of  small  transport  move  and  the  ionospheric  potential  decreases.  As 
time  passes,  these  phenomena  will  begin  to  happen  somewhere  else  in  the 
aurora,  which  thus  shows  a  ceaseless  pattern  of  movement.  It  will  be 
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interesting  to  study  (as  we  intend  to  at  a  later  date)  the  dynamics  with  a 
computer,  to  go  beyond  the  locally  exact  results. 

What  we  have  found  so  far  gives  several  hints  as  to  what  might  happen 
beyond  this  locally  exact  description.  First,  a  slight  generalization  of  the 
hypothesized  spatial  dependence  for  N  and  4>  yields  four  coupled  nonlinear 
ordinary  differential  equations  in  time  which  are  not  obviously  integrable;  if 
they  are  in  fact  not  integrable,  their  solution  (and  a  fortiori  more  general 
solutions  with  no  restrictions  on  spatial  dependence)  will  be  chaotic. 
Second,  linear  perturbation  theory  couples  the  ionization  source  to  the  E  x  B 
-gradient-drift  instability,  and  the  latter,  by  itself,  leads  to  development 
of  spatial  structure  at  finer  and  finer  scales,  through  a  process  of 
bifurcation.  This  might  happen  in  the  full  nonlinear  model  as  well. 

Finally,  a  comment  on  the  role  of  recombination,  which  we  assume  is  small 

2 

for  purposes  of  this  report .  Because  dissociative  recombination  depends  on  N  , 
given  enough  time  it  can  equilibrate  any  ionization  source  which  grows  less 
rapidly  with  N.  However,  this  cuts  the  other  way  too:  in  the  midnight  auroral 
ionosphere  during  the  absence  of  auroral  activity,  N  is  small  enough  so  that 
dissociative  recombination  is  really  quite  small,  and  the  time  scale  is  much 
larger  than  given  in  (1.3).  When  we  speak  of  absence  of  equilibrium,  we  mean 
that  the  time  between  auroral  turn-on  and  recombination  dominance  is  rather 
long  compared  to  the  ionization  time  scale  in  (1.1). 

The  report  is  organized  as  follows:  Section  2  sets  the  stage  and 
introduces  the  model  equations.  Section  3  gives  the  proof  that  there  are 
always  regions  where  transport  cannot  balance  ionization.  Section  4  gives  the 
local  perturbative  analysis,  and  Section  5  gives  the  locally  exact  solutions. 
An  Appendix  gives  a  slight  generalization  of  Section  5's  equations  which  may 
have  chaotic  solutions. 
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2.  PRELIMINARIES 


The  essence  of  the  model  is  to  leave  out  all  phenomena  taking  place  on 
time  scales  much  less  than  tens  of  seconds.  For  certain  effects  this  can  be 
justified,  while  for  others  it  may  be  an  idealization. 

Auroral  electron  (and  even  auroral  proton)  quarter-bounce  times  are 
smaller  than  tens  of  seconds,  so  we  will  not  account  for  delay  times  in 
propagating  particles  from  the  source  region  to  the  ionosphere.  (This  is  not 
true  for  Alfven-wave  bounce  times,  or  wave  periods,  which  are  quite  comparable 
to  our  auroral  time  scales;  it  would  be  interesting  to  study  the  coupling  of 
the  auroral  ionosphere  to  ULF  waves  in  a  model  similar  to  ours.)  We  can  in 
effect  ignore  polarization  currents,  whose  effects  compared  to  Pedersen 
currents  are  in  the  ratio  [Mitchell  et  al . ,  1985;  Cornwall,  1988] 


coc 

-=—  -  0.1-0. 2  (2.1) 
P 

where  0)  =  0.1  sec  1  is  a  typical  frequency,  is  the  height-integrated 

Pedersen  conductivity  (  =  5-10  mho  in  an  aurora)  and  C  *  10  farad  is  the  height- 

M 

integrated  magnetospheric  capacitance: 


C 


M 


(2.2) 


(here  is  the  Alfven  velocity  and  Bj  the  ionospheric  magnetic  field) . 
Moreover,  at  these  low  frequencies  the  magnetosphere  acts  as  an  integrated 
unit  since  wave  and  particle  quarter-bounce  frequencies  are  at  least  as  large 
as  those  that  interest  us.  (Mitchell  et  al.  [1985]  and  Keskinen  et  al .  [1988] 
have  studied  a  case  opoosite  to  ours,  where  no  currents  of  the  type  (1.4)  are 
kept,  but  polarization  currents  are  retained.)  Similarly  we  treat  the 
ionosphere  as  a  unit  in  altitude  by  height-integrating. 

Which  height  region  of  the  ionosphere  we  should  integrate  over  depends  on 
the  energy  of  the  auroral  primaries,  which  determines  their  penetration  depth 
We  will  usually  require 


to  «  v  «  Q 


(2.3) 


where  v  is  the  ion-neutral  collision  frequency,  and  Q  the  ion  (we  assume 


+ 


O  ) 


II 


gyrof requency .  The  first  inequality  allows  us  to  ignore  polarization  currents, 
while  the  second  allows  us  to  drop  Hall  currents  and  to  approximate  the  ion 
drift  velocity  by  that  due  to  E  x  B  drift.  Inequality  (2.3)  implies  that  the 
aurora  does  not  penetrate  below  =150  km  (corresponding  to  primary  energies  of 
only  a  keV  or  so)  ;  for  example,  at  200  km  =  3  sec  \  £2  =  200  sec  .  As 
for  height  integration  of  the  recombination  rate,  this  depends  not  only  on  the 
dissociative-recombination  coefficient  a  but  also  on  the  atomic  ion-molecule 
interchange  rate  at  higher  altitudes,  which  form  the  necessary  molecular  ions 
from  the  preponderant  0+  ions.  Rather  than  introduce  the  extra  complication 
of  ion-molecule  interchange,  we  will  consider  only  two  cases:  predominai.  a  of 
molecular  ions  and  square-law  recombination,  or  linear  and  relatively 
unimportant  recombination  at  higher  altitudes. 

With  these  assumptions,  plus  quasineutrality,  the  equations  to  be  solved 
are  [Cornwall,  1988] : 


V-  (IpV<|>)  -  Q  (<t>  — 4>e ) 


(2.4) 


3n 

dt 


+  ~E 


& 


ro 


2  2 
a(N  -N  ) 

o 


(2.5) 


In  these  equations  all  vectors  are  two-dimensional,  and  we  take  the  x 

coordinate  to  point  north  and  the  y  coordinate  to  point  east;  B  is  constant 

2  2 

and  points  vertically  downward.  The  term  -a(N  -N  )  m  the  continuity 

o 

equation  (1.9)  includes  not  only  square-law  dissociative  recombination  but 

2 

also  any  non-auroral  ionization  sources  (ON  ) .  In  certain  circumstances  we 

2  2  ° 

omit  the  term  -a(N  -N  )  from  (2.5).  In  (2.5),  the  electric  drift  velocity  is 

o 

~E  “  B  (5  X  b>  "  f  dx4,)  (2-6) 

it  is  divergence-free.  As  is  well-known,  the  current-conservation  equation 
(2.4)  defines  a  fundamental  scale  length  A: 

A2  -  I  /Q  (2.7) 


with  A  of  the  order  of  several  ten3  of  km.  Later  we  will  see  that  Vp-ariven 
transport  leads  to  a  new  scale  length  A^,  given  by 
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A 


Nee 


1/2 


(2.8) 


A 


T 


PB  Ip 


which  is  nominally  about  (1/2)  A.  This  is  the  gradient  scale  length  where 
transport  processes  begin  to  compete  with  ionization. 
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3.  ABSENCE  OF  GENERIC  EQUILIBRIUM  AT  a  =  0 


In  Cornwall  [1988]  it  was  shown  that  Eqs .  (2.4,5)  do  have  equilibrium 

solutions  of  a  special  type,  in  which  there  is  a  functional  relation  between  N 

and  (()  (and,  of  course,  is  stationary)  .  Given  such  a  relation,  the  E  x  B 

current  NV  is  identically  conserved,  and  disappears  from  (2.5).  These 
"  E 

equilibria  can  be  two-dimensional  and  depend  on  a  balance  between  ionization 

sources  (first  term  on  the  RHS  of  (2.5))  and  recombination  losses.  Presumably 

other  equilibrium  solutions  exist  in  which  N  is  not  functionally  related  to  <j>, 

but  they  are  not  known  analytically;  in  that  case,  transport  would  also  play  a 

role  in  balancing  sources  and  losses. 

The  situation  is  different  when  recombination  is  unimportant.  We  show 

that  Eqs.  (2.4,5)  with  a  =  0  generically  have  no  time-independent  solutions 

for  an  aurora.  The  reason  is  that  there  must  be  places  where  V  -Vto  vanishes, 

E  ** 

so  there  is  nothing  left  to  balance  the  ionization  created  by  auroral 
primaries . 

An  aurora  generally  has  return-current  boundaries:  two  lines,  running  more 

or  less  east-west,  along  which  <J>  »  $  .  Inside  this  strip — the  aurora  itself — 

we  have  4>  — <J>  >  0.  Outside  the  strip  defined  by  these  lines,  the  ionospheric 

density  N  returns  to  an  ambient  value  N  ,  but  it  is  bigger  than  N  inside.  It 

o  o 

follows  that  for  every  fixed  y,  N  has  at  least  one  maximum  in  x,  where  d^N  = 
0.  This  equation,  9^N(x,y)  =  0,  defines  a  line  in  the  aurora  also  going  more 
or  less  east-we3t.  As  we  follow  this  line  (3ay,  by  increasing  y)  we  must 
encounter  at  least  one  value  of  y  where  9  N  -  0  also,  since  otherwise  N  would 

y 

continually  increase  or  decrease  as  we  go.  So  there  is  at  least  one  point 

where  V  -V  -  0  (because  Vto  «  0),  and  by  continuity  a  finite  region  including 

this  point  where  transport  cannot  balance  auroral  ionization. 


We  can  similarly  show  that  ha3  zeroes  or  near-zeroes.  For  a  generic 
aurora,  V  changes  sign  in  going  from  one  return-current  boundary  to  the 
other  at  any  fixed  y;  this  corresponds  to  the  characteristic  change  of  3ign  of 


E^  seen  in  auroras.  The  equation  V£^(x,y)  «  0  again  defines  a  line  running 
more  or  less  east-west,  characterized  by  an  infinitesimal  tangent  vector  dt . 
Along  this  line  E^  -  0,  and  since  V  x  E  -  0,  we  have 


o 


df-E 


0  , 


15 


(3.1) 


where  the  integral  goes  around  the  auroral  zone.  Since  d£  is  more  or  less  in 
the  y-direction,  (3.1)  implies  there  is  at  least  one  place  where  *  0 . 
Consequently  there  is  one  or  more  places  where  V  s  0  and  transport  is 
negligible  or  zero.  If  recombination  is  small,  again  there  is  nothing  to 
balance  the  buildup  of  ionization  produced  by  auroral  primaries. 

The  equation 

V- (NV  )  =0  (3.2) 

—  ~E 

has,  as  we  have  just  seen,  at  least  one  solution,  say  at  x  =  x  ,  y  =  y  (which 

o  o 

in  general  are  functions  of  t) .  This  equation  is  a  continuous  function  of  x 
and  y,  so  its  zeroes  generically  define  a  line  f(x,y)  »  0  where  transport  is 
negligible.  In  a  finite  region  around  each  point  of  this  line  transport 
cannot  be  in  equilibrium  with  ionization  prodr  ~*-ion,  so  in  general  the  auroral 
region  inside  the  return-current  boundaries  l— s  strips  of  finite  width  where 
equilibrium  is  impossible. 

Given  these  strips,  what  is  the  qualitative  nature  of  the  evolution  in 
time  of  the  basic  equations  (2.4,5)?  We  begin  to  answer  that  question  in  the 
next  section  with  a  study  of  linear  perturbations,  and  then  go  on  to  fully 
nonlinear  and  exact  solutions. 
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4.  LOCAL  PERTURBATION  THEORY 


Local  perturbation  theory  can  be  at  best  a  heuristic  guide,  since  the 
essence  of  auroral  dynamics  is  the  steep  density  gradients  related  to  ioniza¬ 
tion  and  transport  processes.  Nonetheless,  we  will  find  here  a  useful 
introduction  to  the  exact  solutions  of  the  next  section. 

We  linearize  Eqs .  (2.4,5),  writing 


N  =  N  (x)  +  N, 
o  ~  1 


(4.1) 


*  a  ,  .  .  *  i  (k-x-cot) 
0=0  (x)  +  6  e  -  ~ 
o  ~  1 


(4.2) 


The  linearized  equations  are: 


L  +  i  —  [k  0,  d  N  -  k  0,  9  N  +  k  N  d  0  -  k  N  d  0  ) 

1  B  xlyo  ylxo  ylxo  xlyo 


—  0  +  2  ON  N - —  0 

e  l  o  1  e  El 


N 

iN,  k-V0  +  i0  k-^J  -  k^0  N  +  N  0  -  — 0 

1  -  ~  o  l--o  lo  1  o  .2  1 

A 


2  ^El 


(4.3) 


(4.4) 


The  scale  length  A  is  defined  via  (2.7): 

o 


A  -  I  (N=N  ) /Q 
o  P  o 


To  find  the  dispersion  relation,  we  drop  the  inhomogeneous  perturbation  0  . 

El 

Straightforward  algebra  gives: 


(D-k-V  -  -2i0tN  +  i  —  +  -  (k  3  -k  3  )  N 
~  -Eo  o  e  Bxyyxo 

( i  k  •  V0  +  0  ) 

x  - =-=S- - -  °  .  (4.6) 

[-ik-V*I  +  N  (k  +A  )  ) 

~  -  o  o  o 

Here  V„  is  the  electric  drift  velocity  (2.6)  at  0  -  0  . 

-Eo  o 

There  are  several  potential  instabilities  in  this  formula.  First,  for 
long  wavelengths  (kAQ  «  1)  one  finds 
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Im  (0 


(4.7) 


10. 

eN 


.A  V*  ) 
o  o 


2  ON 


The  first  term  on  the  RHS  is  always  positive,  and  is  at  least  as  large  as  T. 


-1 


as  defined  in  (1.1).  These  statements  follow  from  the  equation  (2.4) 

satisfied  by  <b  : 

o 


vV 


O  E 


_1_ 

N 


V<1>  •  VN 
-o-o 


(4.8) 


In  an  aurora  4>  -<b  is  positive,  and  V$  -VN  is  generically  negative,  so 
2o  °  E  -o-o 

A^  v  <j>o  is  positive  and  greater  than  the  potential  drop 

This  ionization-driven  instability  ultimately  will  be  controlled  by 

recombination,  transport,  or  decrease  of  the  parallel  potential  drop. 

Transport  phenomena  are  important  at  large  k,  where  it  can  be  seen  from  (4.6) 

that  the  ionization  term  (~D  is  small.  At  large  k  (kA  »  1)  we  find: 

o 


Ima)-»Y  =  N_1(V  -V-k-V  k •  V)  N 
o  ~Eo  -  -Eo  -  o 


2  ON 


(4 . 9) 


where  k  is  the  unit  vector  in  the  direction  of  k.  The  terms  in  V  can  be  of 

~  Eo 

either  sign,  a  negative  contribution  representing  damping  of  the  ionization- 

driven  instability.  A  positive  term  is  just  a  version  of  the  standard  E  x  B 

gradient -drift  instability  [e.g.,  Linson  and  Workman,  1970].  In  any  case,  the 

magnitude  of  this  sort  of  term  in  Im  0)  is  of  the  order  of  the  inverse  drift- 

transport  time  Td  1  (see  (1.2)).  Simple  geometric  considerations  show  that 

terms  of  both  signs  are  likely  to  be  present  in  a  fully  developed  aurora, 

that  i3,  one  with  significant  E-W  gradients. 

We  can  roughly  estimate  the  cross-over  scale  length  A^  where  drift- 

transport  and  ionization  contributions  to  Im  0)  are  comparable,  by  equating 

these  contributions  from  (4.6)  and  solving  for  k  -  A^  1 .  In  so  doing,  assume 

the  scale  lengths  of  N  and  $  are  the  same,  and  ignore  geometrical  factors. 

o  o 

One  then  finds 


cN  e  i 
o 


Bro 


1/2 


cN  e 
o 


Bf 


Po 


1/2 


(4.10) 


where  in  the  second  equality  we  used  (2.7).  For  the  nominal  parameter  values 
given  in  Section  1,  A^  is  roughly  (1/2) A^,  but  this  depends  somewhat  on  the 
auroral  energy.  Low  energy  auroras  have  smaller  values  for  both  T  and  I 


Po 
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than  high-energy  auroras,  because  they  operate  at  higher  altitudes  and  because 
they  make  fewer  electron-ion  pairs. 

It  is  interesting  to  consider  now  the  inhomogeneous  problem,  in  which  the 
perturbation  0  in  the  driving  magnetospheric  potential  is  retained.  In 

lit  X 

earlier  work  [Chiu  et  al.,  1981]  it  was  observed  that  in  the  static  adiabatic 
auroral  model,  small-scale  (kA  »  1)  fluctuations  in  0  were  shielded  from  the 

lit 

ionosphere,  while  fluctuations  with  kA  >  1  were  not  shielded.  This  makes  it 

difficult  to  understand  small-scale  auroral  ionospheric  fluctuations  as  being 

straightforwardly  driven  from  the  magnetosphere.  Weimer  et  al.  [1985]  have 

confirmed  this  shielding  effect  experimentally.  Making  the  approximations  of 

a  time-stationary  aurora  and  constant  they  used  the  simplified  relation 

which  then  follows  from  (2.4),  or  equivalently  from  (4.4)  by  setting  and 

VM  =0: 

~  o 


<t>Ei(!5> 

0,  OO - T - "  (4.11) 

1+k  A 

o 

(where  <> ^  is  the  Fourier  transform  of  0,  etc.)  and  found  it  to  be  reasonably 
well-satisfied.  However,  at  non-zero  frequencies  this  shielding  effect  should 
diminish.  The  reduction  of  shielding  can  be  seen  even  at  the  low  frequencies 
of  interest  in  the  present  report  ,  by  solving  the  full  inhomogeneous  equations 
(4.3,4).  We  give  the  results  only  for  kA^  >  1: 


01  (k, co) 


0£1 <k,G» 

2  2 
kA 

o 


1  + 


o-k-v  -iy 
-  ~Eo 


r  A  2r°  il 

[iY+-^T!c-V0o]J  , 


(4.12) 


where  Y  (defined  in  (4.9))  is  the  large-k  limit  of  Im  a).  In  the 

approximations  used  by  Weimer  et  al.  [1985],  namely  *  0,  T  »  0,  one  has 

y  *  0  and  so  (4.12)  is  just  the  large-k  limit  of  (4.11)  .  But  when  these 

approximations  are  not  valid,  the  terms  in  square  brackets  come  into  play. 

For  example,  in  an  aurora  with  no  E-W  variation  (k  -  3  »  0)  we  find  that 

y  y 

(4.12)  is,  for  large  k: 


VEo  1 

Vk'(°)  "  ^E! (k'0>)  7^177  — 


which  decreases  only  like  k 


-1 


(4.13) 
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5.  INTEGRABLE  MODEL  OF  TRANSPORT-IONIZATION  COMPETITION 


As  discussed  in  the  Introduction,  when  a  «  0  (no  recombination) 
equilibrium  is  generally  impossible,  because  there  are  regions  where  transport 
is  too  small  to  balance  ionization.  Here  we  illustrate  what  goes  on  in  the 
vicinity  of  a  zero  in  the  transport  term  V^-Vn  in  (2.5)  when  a  =  0,  by 
postulating  that  N  and  4>  in  (2.4,5)  are  no  more  than  quadratic  in  spatial 
variables,  with  time-dependent  coefficients.  Remarkably,  these  nonlinear 
partial  differential  equations  fall  apart  into  four  integrable  ordinary 
differential  equations  in  time,  whose  solutions  show  the  slow  wandering  of 
auroral  features  discussed  in  the  Introduction. 

The  particular  spatial  dependence  used  to  get  these  results  is  not  the 
most  general  quadratic  dependence.  We  briefly  study  the  most  general 
quadratic  case  in  the  Appendix,  and  find  four  nonlinear  first-order 
differential  equations  in  the  time  variable.  These  equations  are  not 
obviously  integrable.  If  they  are  in  fact  not  integrable,  they  may  develop 
chaotically,  an  interesting  possibility  which  will  require  computer  studies. 

Consider  now  Eqs .  (2.4,5)  with  a  »  0.  We  seek  a  solution  of  the  form 


<{>  «=  <|>  +  ^x  +  <|>2y  +  (J^x 


N  =  N  +  N^x  +  N2y  +  N^x 


(5.1) 


(5.2) 


(5.3) 


Here  all  the  functions  with  tildes  depend  on  time  only.  Recall  that  the  x- 

variable  runs  north-south  and  y  runs  east-west.  We  require  that  the  driving 

potential  i  be  positive,  and  that  N  be  negative,  so  that  there  is  a  density 

peak  at  x  *  0  and  a  horizontal  electric  field  with  negative  divergences.  We 

also  require  that  0  >  0  so  that  near  the  origin  (x=y=0)  the  ionospheric 

potential  exceeds  6,  which  has  a  minimum  there.  This  is  a  typical 

£ 

configuration  for  auroral  arcs.  In  a  static  arc  (stabilized,  e.g.,  by 
recombination)  one  usually  looks  for  4>  to  have  a  minimum  and  N  a  maximum  at 
the  same  point  where  6  has  a  minimum;  an  example  is  Fig.  1  (adapted  from 

b 

Cornwall  [1988]).  We  can  arrange  this  as  an  initial  condition  by  setting 
$^(0)  -  N1  (0)  -  0  in  (5.1,2),  but  we  will  see  that  as  time  passes,  the  extrema 
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Potentials  and  Density  at  y  =  0  for  a  Steady-State  Solution  with  Recombination  (after 
Cornwall  (1988]).  The  return-current  region  is  displaced  to  |x|  =  »  ,  where  N  re¬ 
turns  to  its  asymptotic  value  of  0.1  (in  units  of  the  density  at  x  =  0)  and  4>  =  4>r.  The 
distance  scale  is  in  units  of  A„. 


of  (J>  and  N  wander  away  from  each  other  and  from  the  minimum  of  thus 

rearranging  the  relative  strengths  of  ionization  sources  and  transport. 

2 

The  ansatz  (5. 1,2, 3)  can  be  generalized  by  adding  terms  -xy  or  y  ;  thi3 
case  is  briefly  discussed  in  the  Appendix,  but  we  know  of  no  way  of 
integrating  this  general  case. 

Using  the  ansatz  in  Eq.  (2.4)  yields  four  relations.  In  writing  these 
relations  we  temporarily  drop  the  tildes  on  the  variables,  so  that  <{>  and  N 
refer  to  only  the  spatially  constant  parts  of  what  we  usually  mean  by  these 
names.  These  relations  are: 


2<|>3n  +  ni4>i  +  n2<}>2 


(5.4; 


4<t,3N!  +  2Vl  "  “I  *1 

A 

N 

2*3N2  *  “2  *2 
6*3*3  ”  7?  ( W  • 

A 

o 


(5.5) 


(5.6) 


(5.7) 


Here  we  have  introduced  an  arbitrary  density  N  and  length  scale  A  ,  related 

o  o 

by 


A  -  I  (N=N  )/Q  . 
o  P  o 


(5.8) 


Eq.  (2.5)  similarly  yields  four  differential  equations: 


n  +  §  (»2«i  -  V2>  ■  ? 


(5.9) 


"i  *  T  (V3  -  W  '  T  h 


(5.10) 


N  *  —  $ 

2  e  2 


(5.11) 


“3  '  T  ‘VV 


(5.12) 


We  can  eliminate  all  the  (^-variables  in  terms  of  the  N-variables  and  <b  ,  using 

u 

Eqs.  (5.4-7): 
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2A  ,  -  N  x_i  (  A  2N  ,  r  N  \-l  1 

♦  ■  ~r~  *E  1  ■  6Ao  r  )  N  +  ~  +  4Ni  (  ~i  -  2n3  ) 

r\  o'  n  A  * 


4N1  *E 
2A  2N 


o _ 2 

N 


2^-1 

1  -  6A 

O  N 

o 


•r.,  1  ~  6A 

E  l.  O 


N  J 


-1 


(5.13) 


(5.14) 


(5.15) 


(5.16) 


These  equations  can  be  partially  non-dimensionalized  by  agreeing  to  measure 
all  densities  in  units  of  N  and  all  lengths  in  units  of  A  ,  in  which  case  the 

-f  ° 

explicit  factors  of  N  A  in  (5.13-16)  should  be  set  to  1;  N  and  the  N. 

00  1 

(i  -  1,2,3)  are  dimensionless;  and  A  is  the  driving  potential  at  x  =  A  .  We 

£  O 

adopt  these  scalings  in  writing  the  differential  equations  (5.9-12)  in  terms 

of  N  and  the  N. : 

1 


TQ6 


eN 


l  1-6N  J 


TO*  2N 

N  -  - -  - 

2  eN  1-6N 


1-6N, 

2c<b 


3 

4ro$E 

eN 


2(1-2N3, 


(5.17) 


(5.18) 


(5.19) 


N  -  - 


N  N  (1+2N  )  2rQ* 

_ xj  i  c.  o  _ r 

_A  2  (1-2N.) (1-6N  )  +  eN 
BA  3  3  O 

O 


(rk)(“tN22*S:)  (5-20> 


Note  that  N^  is  necessarily  negative,  so  there  are  no  singularities  in 
(5.17-20).  These  equations  are  written  in  the  order  in  which  they  will  be 
solved,  and  are  manifestly  integrable.  First,  (5.17)  is  solved  for  after 
introducing  the  scaled  time  T: 


x  “  I  dt‘  * 

eN  J0  E 


(5.21) 


where  it  will  be  recalled  that  $E<t)  has  the  meaning  of  the  physical  driving 

potential  in  (5.3)  evaluated  at  x  -  A  ,  i.e., 

o 


V*’  '  Ao2  *t(t>  ■ 


(5.22) 


The  solution  of  (5.17)  is  given  in  the  implicit  form 
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1 

6 


N. 


N3  (0) 


n3  +  n3<0)  -  t 


with  (0)  <  0 . 

Combining  (5.17)  and  (5.18)  yields 


N 

N 


2 

2 


N2  =  N2  (0) 


N3(°) 


1/3 


(5.23) 


(5.24) 


Given  N,,  and  N3,  (5.19)  for  can  be  integrated  in  principle;  in  practice, 
this  must  be  done  numerically.  Finally,  with  all  the  N,  known,  (5.20)  is 
integrated  for  N. 

These  last  two  equations  contain  the  squared  ratio 


(5.25) 


of  the  cross-over  length  A  introduced  in  (4.10)  and  the  inverted-V  scale  V  . 

°  -l  ° 

This  is  also  the  ratio  of  the  ionization  rate  T  (see  (1.1))  to  the  nominal 

_1  * 

transport  rate  T^  (see  (1.2)).  We  will  treat  £  as  a  fixed  parameter, 

although  in  principle  it  depends  on  the  auroral  potential  drop  $-6  .  This 

E 

comes  about  because  (see  (2.7))  depends  on  the  Pedersen  conductivity,  which 
depends  on  the  collision  frequency,  which  depends  on  the  average  auroral 
deposition  altitude  in  the  ionosphere,  which  depends  on  the  auroral  energy. 
Since  the  auroral  energy  is  in  part  controlled  by  £,  there  is  an  interesting 
feedback  loop  which  we  can  only  explore  by  going  beyond  the  simplification  of 
height-integration.  The  nominal  value  of  £  is  somewhat  less  than  one,  but  we 
expect  that  £  will  evolve  to  a  value  close  to  one  as  transport  processes  begin 
to  be  important  on  spatial  scales  comparable  to  A^. 

Before  turning  to  the  numerical  solution  of  (5.19,20)  we  note  that  the 
four  differential  equations  can  be  exactly  integrated  in  the  limit  I N3 |  4C  1, 
in  which  case  the  solutions  show  a  linear  instability  of  the  long-wave  type 
discussed  in  Section  4 : 


N3  -  N  (0)  e6T 
N2  -  N2(0)  e2X 

N1  -  N^O)  e4X  +  5n2<0)  (e2T-e4X) 


(5.26) 

(5.27) 

(5.28) 
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N 


=  N { 0 )  e2X  +  N  2 (0)  (l+3$2)  (e4T-l) 

+  \  4n2  (0)  [N  (0)  (0)  ]  (e6T-l)  +|[N1(0)  -  Sn2<0)]2 


(e8X-l) 


(5.29) 


These  solutions  hold  for  arbitrary  and  N2>  but  require  «  1,  and  thus 
N^(0)  1,  X  <  1/6.  Clearly,  this  small-N^  case  is  an  instability  driven  by 
the  auroral  ionization  source. 

Much  more  interesting  is  the  large  IN^I  case,  with  «  -1.  This  case  is 
relevant  when  the  auroral  ionization  source  has  produced  a  central  ionospheric 
density  which  is  large  compared  to  the  unperturbed  density  outside  the  aurora, 
and  there  are  thus  large  gradients  of  density;  they  build  up  the  transport 
terms,  and  allow  the  aurorally  produced  ionization  to  be  carried  away.  We 
study  this  case  by  ignoring  constants  in  terms  like  1-6N2  -»  -6N2,  etc.,  and 
find  for  the  power-law  solutions: 


N 

N 

N 


3 

2 

1 


-X  +  O (tn  x) 


AX1/3  +  0(X'2/3) 


£a  4/3  .  nl  1/3 
2  X  +  0(X  )  , 


(5.30) 

(5.31) 

(5.32) 


where  A  is  a  constant  of  integration,  and  X  now  stands  for  the  scaled  time  of 

(5.21)  plus  the  large  positive  constant  IN  (0)  I  . 

As  for  N,  substituting  the  above  results  in  (5.20)  yields  a  cancellation 

of  the  leading  powers  (which,  if  they  had  not  cancelled,  would  have  yielded 
5/3 

N  -  x  )  .  This  cancellation  is  a  result  of  balance  between  ionization  growth 
and  transport  loss,  which  comes  about  as  the  spatially  dependent  terms  in  the 
density,  namely  the  times  powers  of  x  or  y,  grow. 

Further  analytic  progress  seems  to  be  impossible,  so  we  turn  to  numerical 
calculations.  We  give  results  for  the  following  case: 


N(0)  -  1,  N1(0)  -  0,  N2<0)  -  1,  N  (0)  -  "I,  4  “  1/2  .  (5.33) 

This  corresponds  to  an  auroral  region  where  the  initial  potentials  and  density 
are  symmetric  around  x  ■  0  (recall  that  x  is  a  latitudinal  distance  variable), 
but  have  an  E-W  gradient  as  indicated  toy  N2(0)  and  $2<0)  (■  2/7,  for  our 
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parameters).  We  plot  the  initial  density  and  potentials  at  y  =  t  =  0  in  Fig. 

2.  Here  distances  are  measured  in  units  of  A  ,  density  scales  with  N  ,  the 

o  o 

density  at  x  =  0,  and  the  unit  of  potential  is  the  auroral  driving  potential 

6  evaluated  at  a  distance  A  from  the  origin.  The  return-current  region 
E  o 

(where  4>-<S>  <  0)  sets  in  at  |x|  =  x  =  .882,  at  which  point  N  returns  its 

E  c 

unperturbed  value  of  .22.  The  plots  in  Fig.  2  qualitatively  resemble  both 
real  auroral  profiles  and  the  exact  stationary  solutions  of  Cornwall  [1988], 
as  shown  in  Fig.  1.  (Figs.  1  and  2  cannot  be  qualitatively  compared,  because 
of  differences  in  normalizations,  vertical  and  horizontal  scales,  and  in 
actual  physics.) 

Fig.  3  shows  the  potentials  and  density  at  y  =  0  and  scaled  time  x  =  2  (x 

is  defined  in  (5.21)).  One  sees  that  N  has  increased  by  about  a  factor  of 

two,  and  its  peak  has  shifted  to  the  left,  while  the  ionospheric  potential  <J> 

is  both  decreased  and  flattened,  while  its  minimum  has  moved  to  the  right, 

opposite  to  the  density  shift.  For  other  values  of  y  at  X  -2,  add  1.41y  to  N 

and  0.16y  to  $.  There  is,  of  course,  no  explicit  indication  of  the  generation 

of  shorter  length  scales  in  Fig.  3,  because  the  spatial  dependence  remains 

quadratic  at  all  times.  But  an  indication  of  the  global  breakdown  of  the 

quadratic  ansatz  can  be  seen  in  the  minor  inconsistencies  between  the  position 

of  the  return  current  as  determined  by  =  A  and  by  return  of  N  to  the 

£ 

undisturbed  value  of  0.22.  One  expects,  as  the  perturbative  analysis  of 

Section  4  suggests,  that  generalizing  the  quadratic  ansatz  to  repair  the 

return-current  inconsistencies  will  involve  the  unstable  generation  of  shorter 

length  scales.  However,  a  linearized  analysis  at  X  =  2  is  clearly 

inappropriate;  the  linearized  solutions  in  (5.26-29)  at  X  =  2  are  far  too 

large  to  be  consistent  with  linear  perturbation  theory.  What  we  would  expect 

to  see  in  a  full-scale  numerical  analysis  is  that  major  features  of  the  aurora 

(density  maxima,  potential  minima)  wander  around  through  a  finite  fraction  of 

A  in  times  of  order  T  or  T  ,  or  some  tens  of  seconds,  while  at  the  same  time 
O  ID 

short-wavelength  features  begin  to  grow.  All  this  happens  even  if  the  driving 

potential  <b  remains  constant  in  time,  as  in  our  example.  Naturally,  in  the 
£ 

real  world  4>  also  changes  noticeably  on  these  time  scales,  which  further 
complicates  the  issue.  Our  point  here  is  that  such  variations  in  <J)£  are  not 
necessary  to  produce  auroral  dynamics  resembling  what  we  j.n  fact  see,  and  tnat 
one  cannot  infer  much  about  the  dynamics  of  A  from  the  behavior  of  the 
auroral  ionosphere  on  time  scales  of  tens  of  seconds. 
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Potentials  and  Density  at  t  =  0,  y  =  0  for  a  Time-Dependent  Solution  with  no  Re¬ 
combination  (<}>e  Normalized  to  Unity  at  |x|  =  1).  The  return-current  region  is  at  |x| 
=  0.88,  where  N  has  the  ambient  value  0.22  (in  units  of  ;he  density  at  x  =  0).  The 
distance  scale  is  in  units  of  A0. 


Potentials  and  Density  at  t  =  2,  y  =  0  for  a  Time-Dependent  Solution  with  no  Re¬ 
combination  (4>e  Normalized  to  Unity  at  |x|  «=  1).  The  return-current  region  is  at  |x| 
=  0.88,  where  N  has  the  ambient  value  0.22  (in  units  of  the  density  at  x  =  0).  The 
distance  scale  is  in  units  of  Ao-  The  density  maximum  is  shifted  to  the  left  and  in¬ 
creased,  while  the  potential  is  decreased  and  flattened,  with  its  minimum  moving  to 
the  right. 


6.  DISCUSSION 


In  this  report  we  have  investigated  a  new  class  of  auroral-ionosphere 
instabilities  driven  by  ionization  from  auroral  electrons,  and  involving  a 
parallel  current  which  is  linear  in  the  ionosphere-magnetosphere  potential 
drop.  We  emphasize  that  it  is  not  the  parallel  current  per  se  which  drives 
the  instability,  but  rather  the  ionization  that  it  produces.  The  instability 
grows  on  a  time  scale  of  some  tens  of  seconds  (we  assume  that  this  is  somewhat 
shorter  than  the  nominal  recombination  time  scale,  which  is  true  when  auroral- 
ionosphere  densities  are  not  too  large) .  This  is  comparable  to  transport  time 
scales,  but  we  have  shown  that  there  are  necessarily  regions  in  an  aurora 
where  transport  alone  cannot  relieve  the  buildup  of  ionization,  at  least  on 
the  original  (inverted-V)  spatial  scales  of  the  aurora.  Ultimately, 
dissociative  recombination,  proportional  to  the  square  of  the  density,  will 
stabilize  the  produced  ionization,  but  until  then  the  aurora  responds 
dynamically  by  lowering  the  magnetosphere-ionosphere  potential  drop;  by 
generating  small-scale  fluctuations  which  are  effective  in  transport;  or  by 
moving  the  region  of  low  transport  toward  a  region  of  lower  density. 

A  linear  local  instability  analysis  of  this  problem  is  heuristically 
useful,  but  very  limited  practically  because  of  the  presence  of  zeroth-order 
density  gradients.  We  supplant  the  linear  analysis  with  a  nonlinear  analysis 
based  on  the  hypothesis  of  quadratic  spatial  variations  of  potentials  and 
density.  This  results  in  the  replacement  of  the  original  pair  of  coupled 
nonlinear  partial  differential  equations  by  four  coupled,  nonlinear  ordinary 
differential  equations.  It  is  remarkable  and  unexpected  that  these  equations 
are  completely  integrable.  (A  slight  modification  of  the  spatial  variation 
leads  to  four  coupled  ordinary  differential  equations  which  are  not  obviously 
integrable;  these  may  have  chaotic  solutions.)  The  solutions  to  the 
integrable  equations  reveal  the  expected  transport  phenomena  discussed  above 

In  this  report  we  have  only  discussed  the  case  of  a  time-stationary  driving 
magnetospheric  potential,  which  leads  to  unstable  auroral  fluctuations  on  time 
scales  of  tens  of  seconds.  The  next  step  is  to  consider  the  case  of  time- 
dependent  magnetospheric  fluctuations,  associated  either  with  externally 
driven  ULF  waves  in  the  appropriate  frequency  range  (Pc  4-5;  see  Engbretson  et 
al.  (1986] )  or  as  a  magnetospheric  feedback  response  to  fluctuations  in  the 
auroral  ionosphere. 
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APPENDIX 


Here  we  very  briefly  describe  a  generalization  of  the  integrable  model  of 

Section  5,  which  uses  a  more  general  quadratic  an3atz  for  $,  N,  and  6  rather 

£ 

than  the  specialized  forms  (5. 1,2,3).  The  only  point  to  be  made  is  that  the 
resulting  four  nonlinear  differential  equations  in  time  are  not  obviously 
integrable,  and  so  *-heir  solution  may  be  chaotic. 

The  ansatz  we  use  is 


$  =  <i>(t) 

+  x  .  A .  .  x  . 
i  il  3 

(A.  1) 

N  =  N(t) 

+  x, D  .  ,x  . 
i  13  3 

(A. 2) 

4>e  =  x.$ 

„  .  ,x  ,  . 

Ei3  3 

(A.  3) 

where  A,  D,  and  are  symmetric  2x2  matrices  depending  on  t,  and  the  two- 
E 

dimensional  vector  x^  has  components  =  x,  x2  “  y;  rePeated  indices  are 
summed  over.  Note  that  this  is  not  the  most  general  quadratic  ansatz,  since 
terms  linear  in  x^  could  be  added  to  the  RHS  of  (A. 1,2, 3).  We  need  not  deal 
with  this  additional  complication,  which  only  reinforces  the  lack  of  obvious 
integrability . 

As  in  Section  5,  we  drop  the  tildes  on  the  RHS  of  (A. 1,2, 3)  and  we  scale 

lengths  with  A  and  densities  with  N  (see  5.8)).  Then  in  the  formulas  below 
o  o 

N  and  the  matrix  D  are  dimensionless,  while  A  and  6  have  the  dimensions  of 

L 

potential.  Inserting  the  ansatz  into  the  basic  equations  (2.4,5)  with  a  *  0 
yields : 

Tr  A  -  <|>/2N  (A. 4) 

D  Tr  A  +  {A,  D}  =  j  (A-^)  (A.  5) 

N  -  <}»  (A.  6) 

eN 

o 

D  +  “ (ACD-DEA)  -  (A-6  )  .  (A.  7) 

_ ,  2  N  e  E 

BA  o 

o 

In  these  equations, 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  “architect-engineer”  for  national  security 
projects,  specializing  in  advanced  military  space  systems.  Providing  research  support,  the 
corporation’s  Laboratory  Operations  conducts  experimental  and  theoretical  investigations  that 
focus  on  the  application  of  scientific  and  technical  advances  to  such  systems.  Vital  to  the  success 
of  these  investigations  is  the  technical  staffs  wide-ranging  expertise  and  its  ability  to  stay  current 
with  new  developments.  This  expertise  is  enhanced  by  a  research  program  aimed  at  dealing  with 
the  many  problems  associated  with  rapidly  evolving  space  systems.  Contributing  their  capabilities 
to  the  research  effort  are  these  individual  laboratories: 

Acrophysics  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat  transfer 
and  flight  dynamics;  chemical  and  electric  propulsion,  propellant  chemistry,  chemical 
dynamics,  environmental  chemistry,  trace  detection;  spacecraft  structural  mechanics, 
contamination,  thermal  and  structural  control;  high  temperature  thermomechanics,  gas 
kinetics  and  radiation;  cw  and  pulsed  chemical  and  excimer  laser  development, 
including  chemical  kinetics,  spectroscopy,  optical  resonators,  beam  control,  atmos¬ 
pheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  laboratory:  Atmospheric  chemical  reactions,  atmospheric 
optics,  light  scattering,  state-specific  chemical  reactions  and  radiative  signatures  of 
missile  plumes,  sensor  out-of-field-of-vicw  rejection,  applied  laser  spectroscopy,  laser 
chemistry,  laser  optoelectronics,  solar  cell  physics,  battery  electrochemistry,  space 
vacuum  and  radiation  effects  on  materials,  lubrication  and  surface  phenomena, 
thermionic  emission,  photosensitive  materials  and  detectors,  atomic  frequency  stand¬ 
ards,  and  environmental  chemistry. 

Electronics  Research  laboratory:  Microelectronics,  solid-state  device  physics, 
compound  semiconductors,  radiation  hardening;  electro-optics,  quantum  electronics, 
solid-state  lasers,  optical  propagation  and  communications;  microwave  semiconductor 
devices,  microwavc/millimctcr  wave  measurements,  diagnostics  and  radiometry,  micro¬ 
wave/millimeter  wave  thermionic  devices;  atomic  time  and  frequency  standards; 
antennas,  rf  systems,  electromagnetic  propagation  phenomena,  space  communication 
systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals,  alloys, 
ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  nondestructive 
evaluation,  component  failure  analysis  and  reliability;  fracture  mechanics  and  stress 
corrosion;  analysis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures 
as  well  as  in  space  and  enemy-induced  environments. 

Space  Sciences  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray  physics, 
wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric  and  ionospheric 
physics,  density  and  composition  of  the  upper  atmosphere,  remote  sensing  using 
atmospheric  radiation;  solar  physics,  infrared  astronomy,  infrared  signature  analysis; 
effects  of  solar  activity,  magnetic  storms  and  nuclear  explosions  on  the  earth's 
atmosphere,  ionosphere  and  magnetosphere;  effects  of  electromagnetic  and  particulate 
radiations  on  space  systems;  space  instrumentation. 


